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, We analyse the thermally activated depinning of an elastic string (line tension e) governed by Hall 

dynamics from a columnar defect modelled as a cylindrical potential well of depth Vb for the case of a 
small external force F. An effective ID field Ifamiltonian is derived in order to describe the 2D string 
, motion. At high temperatures the decay rate is proportional to F^^'^T~^^'^ eyip[Fo/ F — U{F)/T], 

' with Fo a constant of order of the critical force and U{F) ~ {eVo)^'''^Vo/F the activation energy. 

The results are applied to vortices pinned by columnar defects in superclean superconductors. 
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I. INTRODUCTION 



o 
o 

u 

^ . 

^ ' The search for new mechanisms of pinning is both scientifically challenging as well as an important problem 
"j^ ' to be studied in view of technological applications of high- Tc superconductors. Even if the transport current is 
. less than the critical one, energy dissipation takes place due to quantum or thermally activated creep of vorticesu. 

^ ' Recent measurements of the critical current and of the magnetization relaxation rate in layered high temperature 
'■^ , superconductors show that columnar defects produced by irradiation with heavy ions can strongly suppress the 

^ vortex motiona. In this paper we present a theoretical study of the thermally activated depinning of a vortex from a 

O '[ columnar defect in the presence of a small transport current (classical creep). 

I ^ I • The inverse life-time F of a metastable state can be written in the form F = Ae~^^'^, where S is the Euclidean action 
I along th&.extremal trajectory and A is the prefactor determined by the fluctuations around the saddle-point solution 
^-H . (see Ref.0 for a general review concerning the decay of metastable states). At low temperatures the saddle-point 
' solution is time-dependent and, consequently, S depends on the dynamics of the system (quantum creep), whereas at 
high temperatures the calculation of S alone does not involve the dynamics (classical creep) . Here, we concentrate on 
high temperatures but go beyond the usual exponential accuracy by calculating the prefactor A, a task which does 
] involve the dynamics of the flux lines as well. 
■ In high- Tc superconductors the dynamics of vortices may be dominated by either the dissipative or Hall term in the 
OO ', equation of motion. Microscopic calculations of the dynamic constanta3 show that the ratio of Hall and dissipative 
0^ ■ coefficients a/77 is approximately equal to uJoTy with loq and r the level spacing between localized Caroli-de Gennes- 
Matricon statesJn the core and the relaxation time, respectively. Recent experimental studies on 90 K crystals of 
YBCO (see Ref.B) have been interpreted as providing evidence for the superclean limit, with wqt ^ 15 below 15 K, in 
which case the contribution of dissipative forces can be neglected in this regime. 

Quantum depinning of flux lines governed by Hall dynamics from a columnar defect has been considered by Sonin 
and Horovitzia for the case oLa small external force. For pancake vortices this problem has been studied by Bulaewskii, 
Larkin, Maley, and VinokurQ. The casCric — J <C jc has been investigated by Chudnovsky, Ferrera, and VilenkinOand 
O ' by Morais-Smith, Caldeira, and BlatterO. On the other hand, the problem of the depinning of a massive string from 
^ a linear object has been solved by Skvortsov in the whole temperature rangeEHI and the thermal depinning of a flux 
• ^ , line governed by dissipative dynamics has been considered by Kramer and Kiiliccj. 

' The present paper is organized as follows: In secti on |l^ we reduce the two-dimensional problem of the string motion 

1 ll^tb " 



^ i to an effective one-dimensional problem. In section |Tjthe decay rate of a trapped string is calculated. In section IV 
- - I the results are applied to vortices in superclean superconductors. 



II. ID EFFECTIVE HAMILTONIAN 



Let us consider a string which is pinned by a columnar defect in the presence of an external force. Both, the 
cylindrical defect as well as the vortex are directed along the c-axis of the anisotropic superconductor and the external 
magnetic field is supposed to be sufficiently small, such that the interaction between the vortices can be neglected. 
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Furthermore, we consider the situation where each vortex is pinned by an individual defect, i.e., the concentration of 
defects is assumed to be larger than that of vortices. The free energy density of the string describing this situation is 



e / du 



2 



G[u]-2(^) +U{u), (1) 



where 



Uiu) = V,yi [^ul + ul ) + Ve^tM. (2) 

Here, e is the elasticity of the string and and Uy describe its displacement along the x and y directions, with 
the z-axis chosen parallel to the defect. In Eq. (||), Vcyi {\/ux'^ + Wy^) denotes the cylindrical pinning potential and 
Vext = —Fux is the forcing potential. The function Vcyi (r) is supposed to be monotonously increasing and restricted 
from below and above. We assume boundary conditions u (±L/2) — 0. 

If F 7^ 0, the state of the vortex becomes metastable. Our main goal is to investigate the decay rate as a function 
of temperature T and force F, where the external force is assumed to be small. The Lagrangian of the vortex (in real 
time dynamics) can be written as 

/•+i/2 / \ 
L[ux,Uy\ = J dz l^auy-^ - G[ux,Uy]j , (3) 

where a is the Hall coefficient. The corresponding equations of motion take the form 

dux _ dU_ _ d'^Uy 

dt ~ duy~^ dz^ ' ^ ' 

duy dU d'^Ux 

Eqs. d) and §) can be formulated as the equations of motion of the ID Hamiltonian density 



H = 



U 



(^'a)^2(9z) 2a^{dl) 



(6) 



where we have used the definitions x = Ux and p = auy. Using the variational procedure for the Hamiltonian density 
H we obtain 



dH _ d_dH_ _ dU{x,pla) _ e d'^p 
dp dz dpz dp dz^ ' 



(7) 



dU d dH dU(x,p/a) d^x 

'^-^'-d-zd^^r-^df^^'di^- 

and one can easily see that Eqs. (|^) and (j^) are equivalent to Eqs. (^) and (j^) with x — Ux and p — aUy. Thus we 
have reduced the 2D problem of the motion of a vortex governed by Hall dynamics with an action given by (|3|) to 
the ID problem of the motion of a vortex descrUjied by the effective action J dzdt {px — H) , with H given by Eq. (||). 
The above idea has been introduced by VolovikO for the tunneling of 2D vortices in a liquid helium film, the motion 
of which is equivalent to that of a_massless particle in a magnetic field. Later, this idea has been used by Feigel'man, 
Geshkenbein, Larkin, and LevitEj within the context of vortex tunneling in superclean high-Tc superconductors. 
Chudnovsky, Ferrera, and VilenkinEi have generalized this idea for the description of the depinning of a flux line 
governed by Hall dynamics from a columnar defect near criticality. In their case U{ux, Uy) = — + cUy/2 

and the elasticity along the y-direction can be neglected if jc — j ^ jc, i- e., the problem can be reduced to the ID 
massive string problem. Above we have generalized the problem for the case of an arbitrary potential and nonzero 
elasticity in the y-direction. 
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III. DECAY RATE 



The decay rate F of an arbitrary metastable Hamiltonian system at high temperatures is given by the expressionO 

huiQ ImZ 

= — , 9 

TT Z 

where luq is the unstable mode growth rate and Z is the partition function of the system under consideration. Eq. (|^) 
is apphcable if the temperature T satisfies the condition T ^ hujo^ otherwise quantum effects become relevant. If the 
transition from quantum to classical behavior in the decay of the metastable state is of second order, the parameter 
To = hujQ/2'K determines the temperature of the crossover. Howeverpin general the transition can be of first order, in 
which case the crossover temperature differs from hLJo/2TT (see Ref.EJ). Note that at low temperatures the equation 
for the decay rate can be written as HT — 2TlmZ/Z. We emphasize that Eq. (|^) can be applied only if the system is 
properly equilibrized, i. e., if its characteristic relaxation time is much smaller than T~^. 

We assume that the quasiclassical approximation is applicable. With the partition function written as a path 
integral, 



Z ^ {Vu} exp I I , 



(10) 



we can use the steepest descent method for its calculation. In this approximation the partition function is determined 
by its stationary points. The most significant contribution to Z arises from the trajectory u(z, t) = — the position of 
the minimum of the potential, whereas the imaginary part of the partition function is determined by the neighborhood 
of the saddle-point solution u(z,t) = Uo(2;) which is time-independent at high temperatures. Hence, the problem is 
reduced to the calculation of the unstable mode growth rate ujq , and of the real and imaginary parts of the statistical 
sum. Below we shall calculate all these quantities. 

A. Unstable mode growth rate 

At high temperatures the saddle-point solution does not depend on time. The stationary extremal trajectory Uq{z) 
for the Euclidean action 



r+L/2 


l-+h/2T , 


1 dz 




I-L/2 J 


-?i/2T \ 



' dr 



(11) 



corresponding to the real-time Lagrangian (||) satisfies the Euler equation 



az^ aw. 



The unstable mode growth rate then is determined by the negative eigenvalue of the operator (5u5'euci|uo(z)- Near the 
saddle-point solution uo(2:) the perturbed solution can be expanded in the form 

Uxiz, t) = u^oiz) + i>{z) exp {i2TrTT/h), (13) 
Uy{z,T) — ip{z)e:sq){i2TTTT/h), (14) 



with small distortion amplitudes 4'{z) and (p{z) for T < Tq. Substituting (^3|) and ( p^ into the equations of motion for 
the Euclidean action ( |ll[ ) and expanding, we obtain the system of equations determining the crossover temperature 
To 



d'lp , , d'Vcyl 



'Iz^ ^ ''''''' • dr^ 



ctz^ r dr 



= + aujo'^ = 0, (15) 

r=Ua;o(z) 



= Hyip — aajQip = 0, (16) 
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subject to the boundary conditions ip{±L/2)^ (/?(±L/2) = 0. The (lowest) eigenvalue luq is related to the crossover 
temperature Tq through the condition ljq = 211X0/ h. 

The operator has only one negative eigenvalue: it can be easily seen that in the limit L ^ 00 the function 
duxo/dz is an eigenfunction of the operator with zero eigenvalue. Its derivative du^o/dz can be understood as 
the "velocity of a particle" moving in the potential U{uxo, Wyo) with the "velocity" changing its sign once. Hence, 
duxo / dz is an eigenfunction of the ID Schrodinger operator with one node and there must be another function 
corresponding to the ground state "wave function" with a negative eigenvalue. We conclude that the operator has 
one negative and one zero eigenvalue. The operator Hy is positive and has only positive eigenvalues. The "potentials" 
for both the operators Hx and Hy are plotted in Figs. 1 and 2. 
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Fig.l The "potential" Ux(z) for the ID Schrodinger operator Hx (see Eq. (|l5|)). The characteristic size of the potential is given by 
D = \/2eVb/-F and its height is equal to re. 




Fig. 2 The "potential" Uy{z) for the ID Schrodinger operator Hy (see Eq. (^^). The characteristic size of the potential is given by 
D = ^/2eVo/F and its height is equal to re. 

Eliminating the function ip from the system of Eqs. (|l5|) and ( |l6| ) we obtain (we denote eigenvalues of HxHy by A 
and those of Hx by A) 



HyHxii^ = AV'. 



(17) 



The unstable mode growth rate loq is determined by the negative eigenvalue A-i of Eq. (|17l), A_i = — (awo) • In 
Appendix ^ we derive a variational principle for the present problem and show that the negative eigenvalue of 
Eq. ( p^ indeed exists and is unique. (We exploit that the operator Hy^^ HxHy^^ is Hermitian and has the same 
eigenvalues as HyHx). Furthermore, we derive both upper and lower bounds on the eigenvalue A_i of the operator 
HxHy in Appendix such that we finally arrive at the estimate 



A-il 



2'Ka 



2^~aVo'' 



(18) 



with <^ e [1.1530,2.7314] a universal constant and Vq denotes the depth of the potential Vcyi. The result ( [l8|) holds 
independently of the details of the pinning potential Vcyi as long as the driving force F is small, F <^ F^- We note 
that Tq ^ F'^; the same depencLeace holds for a string with dissipative dynamics. For a massive string depinning from 
a linear defect Tq ~ F, see Ref.El. 
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B. Ratio ImZ/Z 



The most significant contribution to ReZ arises from the neighborhood of the minimum of the potential Vcyi ■ In 
this region we can write (see Eq. (||)) 

t/(x,p/a)-|x2 + ^p^ (19) 
where k = (PVcyi/du'^\u=o- The equations of motion (0) and (^) take the form 

a^x^ L~e^)p=HoP, (20) 



-HqX. 



(21) 



Hence, p = ci^Hq and the Euchdean action of the vortex near the equihbrium position can be written in the form 
(we make use of Eq. (0)) 



Si 



1 



Eucl 



j-+h/2T 


r+L/2 , 


/ I 


dz X 


l-h/2T J 


-L/2 \ 



(22) 



The same procedure as above provides the variation of the Euchdean action near the thermal saddle-point solution 

X, (23) 



'-'Eucl — 2 



r+h/2T 


r+L/2 , 


/ dr 


' dzSx { 


'-h/2T J 


-L/2 V 



(92 „ 



where the operators and Hy are given by Eqa_(|T^) and (|T^). Using 
in the path integrals for S^^^^ and S^^^^ (see Ref.li3), 



the usual measures for harmonic oscillators 



and 



dfi^ = 



dot ( -^'g^H^' 



1/2 



n 



V2Trh 



1/2 



2Trh 



(24) 



(25) 



(here and are expansion coefficients over the systems of normalized functions) and defining the dimensionless 
operators Hq, h^, and hy as Hq/k^ Hx/ k, and Hy/ n, respectively {k, — V/yi (0)), we obtain for hnZ/Z 



ImZ L 



2 V 2ttT 



l^xy,o{L) 



L/2 
1/2 



1/2 

oz J 



1/2 



(26) 



where (L — > oo) 



+L/2 



U = e 



(27) 



-L/2 



is the activation barrier and L is the length of the string (we performed standard integration over the shift mode, 
see Refs.ll3lll). The prime in Eq. (26) indicates that we exclude the zero eigenvalue of the operator —a^Hy^-^ + 
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Hx- fJ'xfiiL) and fJ-xy.oiL) are the "zero" eigenvalues of the operators hx and hy^^hxhy^'^ defined on the interval 
[-i/2, +L/2], L ^ oo. 



The eigenvalues of the operators —^-^ + and — + hy^hxhy^ are equal to 



/27raT 



Mo,a and 



/27raT 



f^ + fJ-xy,b, n = 0, ±1,±2. 



(28) 



with fiQ^a and the eigenvalues of the operators Hq and hl^'^hxhy^'^ , respectively. Substituting these eigenvalues 

into Eq. ( |26| ) and calculating the product over n using Jl^i (l + /v?^ = sinh {ttx)/ttx we obtain 



ImZ aL 



Z 



T 



Ah V 2tik 



exp 







/ ' dz\ 


J-L/2 ^ 





1/2 



u 

T 
exp 



1/2 



Hk 
2^ 



^^xy,0 {L)_ 



a 

n" [i-cxp(-ifvTO) 

b 



(29) 



The double-prime sign reminds us that we excluded the negative and zero eigenvalues of the operator hy hxhy . 



C. Result for V 



"1/2" "1/2 

The positive part of the spectrum of hj hyhj consists of M discrete eigenvalues iixyfi and a continuous part 
l^xy(k) = (l + tk'^ / kY with a spectral density p^{k). The spectrum of is continuous, /io(fc) = (l + ek'^ j k)^ with 
the spectral density p'^{k). Using 



n/Wa) 

a 

n/(Mxy,6) 

6 



exp 



dk 5p{k)\nf {p{k)) 



with 5p(k) = /o'^(fc) — P^ (k) we arrive at the final expression for the decay rate. 

U* 



r = 



aTn / T 



2n^ V 27rK 

/ OO 

exp i J dk Sp{k) In 



expl-- 





/ ^ dz\ 

J-L/2 ^ 


'duo' 


i 

) 


1/2 


^2^,0 (i) 


/ sin^ 








.Pxy,Q [L] _ 



— (-^(-^'^'))Jljn 



M r 



1 - exp ( -^y/Pxy,b, 



with 



U* = U 



2a 



I dk sp{k) (i + -^fc^) - 



(30) 



(31) 



(32) 



the quantum renormalized activation energy. Let us investigate the correction to the activation energy ( |32[ ) arising 
from large fc-values. We will see that the large k modes lead to an ultraviolet divergence which we have to cut off by 



some physical length scale. It is convenient to rewrite the integral / dk' Sp{k') {l + ek'^ / as 



No{k) 



Ns{k) 



j dk'5p{k'){v + ^k'^)^ J2 {l + ekl/n^)- J2 (l + ^^n^')' 



(33) 
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where fc„ and qn are trivially related to the eigenvalues of the operators Hq and h\/'^hxh]/'^ and No{k) and Ns{k) 
denote the number of eigenvalues inside the interval [0, k]. Note that Nq{oo) — Ns{oo) = M + 2. The two operators 

/ig and hkj'^hjh^j'^ define two scattering problems. We define appropriate scattering states ipkiz) and i'qiz) with 
asymptotics e*'^^ and e*''^ {z —oo) and g^''^ ^ g^i^+^ii) (^z +00), with S{q) the appropriate scattering phase 
shift. The general solutions of the "Schrodinger" equations can be written as "ifkiz) = Ciipk{z) + C2ip-k{z) and 
'i'q{z) = C[ipq{z) + C'2^-q{z). Since we require that the "wave functions" satisfy jtbe condition ^p{±L/2) — 0, 
L ^ -\/eVo/ we obtain the discrete levels through the phase equations (see also Ref Jl3) 



knL = Trn, q-nL + S{q) = irn, (34) 



and hence 



,0 o /7rn\2 /tttt, — (5(q)\^ 27rn , , 2q ^, , , , 

fc' -'?«-( — ) -[ ^ +j^S{q) :^ +^5{q). (35) 

For q ^ \/K/e we can make use of the semiclassical approximation for the operator IhJ'^hJritl'^ and calculate the 



phase shift directly: With the usual Ansatz '^{z) ^ e I '^^^^ '^^ the function q(z) can be found from the equation 

{eq^{z) + {eq\z) + Uy) = K^Eq, (36) 

where i?g = (l + tq^ j nf' ~ e^q^ / is the eigenvalue of the continuous spectrum of h^ij'^lrij'^hkj'^ ■ At large 5 we obtain 




g(z)^,/-ii;V^|l- -- |. (37) 

Similarly, for the operator /iq we obtain 

k(^z)^^E]!'\\-^\^ (38) 
and we arrive at the difference in the phase shift 



-L/2 

Taking into account that at large g, ~ ^q^ j and that {/^^(z), C/j^(^;) <C k for |z| < v^2eVo/F, we obtain 



(40) 



i.e., k"^ — q^ — 2n\j2V{^l tj FL. At large k the number of states in the interval dk is equal to Ldk/2TT and the integral 



fe 



,2\ JVo(fc) 



1 + U 1 ^ (fc^ - ,^) ^ I dfc' (41) 

diveijges-Jinearly at large fc (for the case of a massive string the correction to the activation diverges as Infc, see 
Refs£3'E3). The integral then has to be cut off at some wave vector k* : For a string the natural cutoff is 7r/r, with r 
the radius of the string. For a vortex parallel to the c-axis of an anisotropic superconductor k* ~ 7r/max (^cj d)^ with 

the coherence length in the c-direction and d the distance between the superconducting layers. 

Cutting off the integral in (^) at large momentum k* we obtain 

The theory is self-consistent if the correction to the activation energy is small as compared to U (here we have assumed 
the cutoff k* to be large such that we can neglect the contribution from the bound states). 
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In case of a continuous transition from quantum to classical behavior the result Eq. (E^) with Eq. ( ^2| ) is applicable 
for any T > Tq except for a narrow temperature interval ^ h^^'^ around Tq (see RefS). An explicit expression for 
the decay rate can be obtained only for the case ?i — > 0, i. e., in the classical limit. In this limit we can expand the 
exponents in (|3^), 1 — exp{—'hK/aT) ~ Tin/aT and the quantum correction to the activation energy (|3^) tends to 
zero. Performing this steps we obtain the result 



2'Ka 



+L/2 



dz 



L/2 



dz J 



1/2, 



-,1/2 



det 



det' 



1/2 



(43) 



Using the same technique as before, see Eq. (Eq), we arrive at the more suitable form 



kL 

2'Ka 



2ttT 



exp 



det 



-L/2 \dz J 

1/2 



1/2 



det 



det ( hf) 



det 



{hy) 



1/2 



(44) 



where pLxfi{L) is the zero eigenvalue of the operator hx{L 
The caJculation of the instapton determinant ratios ([ 



can be elegantly performed using the Gelfand-Yaglom 



formulaO. The renormalizedu determinant delH of the Schrodinger operator H = —d^/dz^ + p{z) is the solution 
of the differential equation Hf{z) = with the initial conditions f\z=-L/2 = and dfld?]- r/n = 1. The value of 
the function f{z — L/2) provides the renormalized determinant. As has been shown in Ref.Ej, in the limit ^ 0, 
/ix,o(i) (Fa/yo)exp - ./i^e [L - 2^/2^ / F) , det(/io)/|det(/i:r(i))| exp^i/^L^, and det(/io)/det(/ij,) - 
(Fa/Vo) exp {2^2kVo/ F) , where the coefficients of proportionality depend on the detailed form of the pinning potential 
and a is the characteristic radius of the pinning well. The eigenvalue fixy.-i = ^^'^F^/i^'^Vo, see Eq. (p^), and we 
make use of Eq. ( ^ . Substituting these dependences into Eq. ( ^ we obtain 



VoL 



5/2 



1/2 



exp 



U 
T 



F 



(45) 



where 3 is a dimensionless proportionality coefficient which depends on the detailed form of the pinning potential, 
T = ay/eVo, Fc is the critical force, and the activation energy U is given by Eq. (p7|). As an illustration we have 
explicitly solved the differential equations for a pinning potential given by the equation 



cyl 



1 



cos^ — 
\2a 



I) 



(46) 



where 6{t) is the step-function (9{t) 
Fei\u,\-a). 

We find fi^^oiL) = (24i^a/7ryo) exp 



1 if t > and 9{t) — otherwise). The external force has the form F(u^) 

det(/io)/|det(/i:r)| = (1/2) exp 



and 



/K/e [L - 2^/2eVQ/F) 

det(^o)/det(/ij^) ~ 4Fa/7rVb exp (2v^2kV^/F), and substituting these expressions into Eq. ( |44| ) we obtain g — 
7r^/25/4 ^ 1 32^. 

A real physical system will involve a finite cutoff k* , which we can account for in Eq. ( |3l| ) by a restricted integration 
over the modes k < k*. However, in making use of the Gelfand-Yaglom formula we actually account for unphysical 
fluctuations with k > k*. In order for the classical result ( psf ) to be valid we then must require the cutoff k* to be 
sufRciently large such as to validate the Gelfand-Yaglom approach. On the other hand, k* has to be small enough in 
order to arrive at a small quantum correction in Eq. (135 



IV. CONCLUSION 



Let us first specify the regime of applicability of our results. The zero temperature Euclidean action takes the value 
•SeucI {T = 0) = ail, where O is the volume encircled by the string in the course of the tunneling motion. In the x- and 
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y-directions the characteristic size of the loop is Vq /F and the length of the string segment along the ^-direction is of 
order y/eVo/F, hence 5euc1 (T = 0) - a^/7JVo{Vo/ Ff . On the other hand, 5euc1 {T = Tq) = HU/To ^eucI (T = 0) , 
i.e., even if the transition from quantum to classical behavior at Tc > Tq is first order like, its temperature will still 
be of order ~ Tq and our results are applicable in the regime T > Tq. 

Next, we estimate the temperature Tq using parameters for the moderately anisotropic superonductor 
YBa2Cu307_5. We choose our columnar defect in the form of a cylindrical cavity of radius a ~ Cafc(O), wherej^pf, is 
the coherence length in the a6-plane. The expression for the pinning potential then can be written in the formtj 



whera-so = ($o/47rAa6)^ is the relevant energy scale and the factor v {0 < v < 1) describes the pinning efficiency 
factoinJ, $0 and Xab are the flux quantum and the penetration length in the afo-plane, respectively. For the model 
potential (^ ) the temperature Tq, the critical force Fc, and the temperature hn/a can be easily calculated (the 
Hall coefficient a is related to the electron density n as a = T:hn at T = in the superclean limit luqt 3> 1): 
To - (!^C/647r2) (eo/n^l^) {F/F,f, F, = {V^iy/U) eo/Cab, hi^/a = lyheQ/AaCb- Substituting U = 16A, Xab = WOO A, 
Xc/Xab = 5, a = TThn, n = 2-10^^ cm~'^, e = cqX^i^/ Xl^and ^ ~ 2, we obtain Tq = 0.6K- ^{F/Fc)^, and Tin/a — 15K-iy. 
Takiug into account that ly might be ~ 0.1, see RefB3, we see that at T ~10- 15 K when the Hall force appears to be 
largeH, the conditions T '3> Tq, Tin/ a are well jSatisfied. _ 

The ratio a/rj 15 for T <i 15K (see Ref.|j) was obtained from indirect measurements. In Ref.Eil direct measure- 
ments of the Hall angle in 60 K YBa2Cu307_a are reported to yield a ratio a/77 ~ 1. We wish to point out that even if 



a/?7 1 it is still possible to use Eq. (45): We have shown that the inverse lifetime of a vortex pinned by a columnar 
defect behaves at high temperatures as 

r^i^^/^T-V.exp(-^ + ^). (48) 

It is interesting tp-,note that the same dependence of the decay rate on the external force and temperature has been 
obtamed in RefsO for a flux line governed by dissipative dynamics. Moreover, the equation for the decay rate from 
RefsH can be obtained up to a numerical prefactor from Eq. ( ^ ) by the substitution a ^ t], indicating that in the 
regimes a < 77 or a ^ 77 it is possible to use Eq. ( ^5| ) with a — > max (a, 7;) . 
The scaling law (Es) for the decay rate leads then to the resistivity scaling 



piF) 



|i^exp(2v/2^/F 



exp(-^V (49) 



where the factor in square brackets represents the correction to the standard result, which arises from the inclusion 
of classical fluctuations around the saddle point (prefactor). 

We note that Eq. (^) is not valid for T > FU /2\/2kVq as the quasiclassical approximation is not applicable at 
these large temperatures. It is well-known that the problem of a vortex pinned by a cylindrical potential is equivalent 
to that of a 2D quantum particle trapped in a radially symmetric 2D potential. The latter always has a bound state 
in 2D, and thus the vortex is pinned by the defect at any temperature. However, the tiiewjial fluctuations lead to a 
large downward renormalization of the pinning energy at sufficiently high temperaturef3l23. 

Finally, let us make an estimate of the cutoff momentum fc* . The theory developed above works well if fc* ^ k/e; 
with fc* = Tr/d, where d — 12 A is the distance between superconducting layers, we obtain fc*/-\/ k/c ~ 5, such that 
this condition is marginally satisfied. 

Briefly summarizing, we have considered the problem of the thermally activated depinning of a flux line governed 
by Hall dynamics from a columnar defect in the presence of a small (F <C Fc) external force. We have shown how 
to reduce the 2D problem of the string motion to a ID effective problem. The expression for the decay rate has 
been obtained for the whole temperature region where the saddle-point solution is time-independent (see Eq. ( |3l| ) 
and Eqs. (Q) and (^5|) for high temperatures). An analytical expression for the classical asymptotics of F has been 
calculated for a model potential as given by Eq. (^6|), and we have discussed possible applications of these results to 
high-Tc superconductors. 
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APPENDIX A: EXISTENCE AND UNIQUENESS 



We prove that a solution of Eq. with a negative eigenvalue exists. Second, we show the uniqueness of this 
solution. We mention that in the limit L oo the boundary condition ■0^0 for |z| — s- oo is asymptotically satisfied. 

Existence. As we already know, the operator has one negative eigenvalue. We show that the operator HyHx 
has a negative eigenvalue, too. It is possible to rewrite the equation HyH^ip = Xip in the form 

hIh,hI{hP^) = \{HPi^) (Al) 

(the operator Hy has only positive eigenvalues, hence, operators like Hy ^ or Hy ^ are well defined and Hermitian). 
The operator Hy H^Hy is also Hermitian and its lowest eigenvalue can be obtained by minimization of the functional 

(xlx) 

in the space of functions with integrable modulus squared. Using the fact that Hy'^ is self-conjugate and defining 
1 

(p = Hy ^ Xi tti6 functional can be rewritten as 

F[<t>] = - . (A3) 

{Hy-(ly\Hy--4>) (c^\Hy | 0) 

The operator H^ has a negative eigenvalue, i. e., there is a function \(j>) such that {(l)\Hx\<j>) < 0. The operator Hy has 
only positive eigenvalues, i. e., the form {4>\Hy \4>) is always positive. Consequently we constructed the Rayleigh- 
Ritz principle for the equation HyHx^ = Xip and proved that there is a function with a negative average value. The 
eigenvalue A_i is even lower, thus we have demonstrated the existence of a negative eigenvalue for the problem (p^). 
Uniqueness. Let us rewrite the operator H^ in the following form 

i?. -^A„|a)(a|- ^ A«|a)(a|+A_i|-l)(-l|, (A4) 

a Q^ — 1 

where Aa are the eigenvalues of Hx, A_i is the lowest eigenvalue of the operator Hx, the vector | — 1) denoting the 
"ground state" of the operator Hx- The operator Hy HxHy can be rewritten in the form 

hIHxHI = hI X^\a){a\Hl +hI\-i\ - 1){-1\hI =A + B. (A5) 

The image of the operator B = Hy X-i \ — l){~l\Hy is one-dimensional, i. e., S is a self-adjoint projector of rank 1. 
The operator 

^ = ^1 Aa|«)(a|-l-0-|-l>(-l|j hI (A6) 

is a nonnegative Hermitian operator with two zero eigenvalues. We perturb the operator A with a projector of rank 
1 and make use the following theoremc^ 

Theorem. Let A be a self-adjoint operator and B a self-adjoint operator of rank 1. Take an interval D on the real 
axis and denote by n{D) the number of eigenvalues of A in this interval. Then the number m{D) of eigenvalues of 
A-\- B within D satisfies the estimate: 

n{D) - 1 < ni{D) < n{D) + 1. (A7) 

Applying this theorem to the interval (— oo,0) (we exclude zero!) one can easily see that the operator Hy HxHy has 
not more than one negative eigenvalue. But the eigenvalues of this operator are the same as those of the operator 
HyHx so that our problem has only one negative eigenvalue. 
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APPENDIX B: NEGATIVE EIGENVALUE PROBLEM 



We construct upper and lower bounds for the negative eigenvalue A_i of the operator HxHy. In the limit _F ^ the 
lowest eigenvalues of the operators Hx and Hy do not depend on the detailed form of the pinning potentiaS. This is 
due to the rapid decay of the eigenfunctions in the region \z\ > D, where the form of the potential is irrelevanL, The 
lowest eigenvalue of the operator is equal to A_i = —fi^F'^ /2Vo, with fi the root of /ztanh/x = 1 (see RefsO). In 
the region jzj < I? the normalized "ground-state" wave- function of the operator is 



cosh 



0-1 w 



D cosh 



D 



(Bl) 



As the characteristic depth k of the "potential well" Uy{z) is much larger than the "size quantization energy" e/D^, 
see Fig. 2, we can approximate the eigenvalues and eigenfunctions of the operator Hy by those of the infinitely deep 
quantum well 



, 1 fnn vr 

0„(.) = -=cos(— z + 



1 



1,2. 



Using the variational principle (A3) we can write 



< 



2Vo 



(B2) 



(B3) 



Substituting the eigenvalues and eigenfunctions (B2) into Eq. (B3) we obtain 



< 



F4 

= -1.3292-^. 



(A<2+7r2(2m+l)V4) 



In order to obtain a lower estimate A_i > B > ~oo, (see (A3)) we use the following inequality 

A-i 



where the functions belonging to V satisfy the conditions i) 



A_i = min0e£2F[0] > min^ep , , 

{(l)\Hy \q 

-- 1, m) 



(B4) 



(B5) 
z). The 



< 0, and ill) (t>{z) = 

last condition is a consequence of the parity symmetry of the operator H^Hy. Following the discussion in section [II A, 
the operator has one negative eigenvalue with an even eigenfunction. For any odd function (p, {(l)\Ha:\(l)) > 0, thus 
an odd function cannot minimize the functional (A3), hence the eigenfunction corresponding to the lowest eigenvalue 
of the operator HxHy has to be even. Our goal is to find some (large) constant b such that {'^\H~^\'^) > 6 > for 
any 'i> £ V. We rewrite 



\H- 



E 



> \{m)\ 



(B6) 



and expanding (j) in the eigenfunctions of the operator Hx , — 



c_i0_i -1- (p, we can write 



As {(l)\Hx\(f>) < 0, — c^]^|A_i| + cfXi + C2A2 + . . . < 0, and taking into account that = 1 we obtain 



(B7) 



1 > |c_i| > 



' Ai 
lA-il + A] 



= 0.9191, 



(B8) 
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where we used the fact that A_i = -1.4392e/D2 ^ 7.8309e/D^ (see Refs.li°IO). As \{4'\4'[)\ < \\4>\\ \\4>[\\ < 

1 - c2_^ = 0.3941, we obtain > f{\c^i\,K), where 

fi\c-i\,K) = cl,\{cj^^M\' + K'~ 2\c^,\\{cj^^M\K. (B9) 

We note that |c_i| G [0.9191, 1] and K e [0, 0.3941]. A simple analysis shows that the minimal value of /(|c_i|, K) is 
equal to /(0. 9191, 0.3941) 0.119 and substituting this value into Eqs. (E^) and (B6) we obtain 



4 



A-i > -7.4603^. (BIO) 



Vr 



Hence, we have shown that in the limit F — > the lowest eigenvalue A_i of the operator HxHy, A_i behaves as 
~jF^/Vq with 7 e [1.3292,7.4603]. The temperature Tq and the eigenvalue A_i are connected by the equation 



(Bll) 



Tq = hy \X^i\/2'Ka and we obtain the estimate 

° 2n aVo 

with C e [1.1530,2.7314]. 
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